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ABSTRACT 

The Lyman-a emission line dominates the far-ultraviolet spectra of late-type 
stars and is a major source for photodissociation of important molecules including 
H2O, CH4, and CO2 in exoplanet atmospheres. The incident flux in this line 
illuminating an exoplanet 's atmosphere cannot be measured directly as neutral 
hydrogen in the interstellar medium (ISM) attenuates most of the flux reaching 
the Earth. Reconstruction of the intrinsic Lyman-a line has been accomplished 
for a limited number of nearby stars, but is not feasible for distant or faint host 
stars. We identify correlations connecting the intrinsic Lyman-a flux with the 
flux in other emission lines formed in the stellar chromosphere, and find that 
these correlations depend only gradually on the flux in the other lines. These 
correlations, which are based on HST spectra, reconstructed Lyman-a line fluxes, 
and irradiance spectra of the quiet and active Sun, are required for photochemical 
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models of exoplanet atmospheres when intrinsic Lyman-a fluxes are not available. 
We find a tight correlation of the intrinsic Lyman-a flux with stellar X-ray flux 
for F5 V to K5 V stars, but much larger dispersion for M stars. We also show 
that knowledge of the stellar effective temperature and rotation rate can provide 
reasonably accurate estimates of the Lyman-a flux for G and K stars, and less 
accurate estimates for cooler stars. 

Subject headings: exoplanet: atmospheres — stars: chromospheres — ultraviolet: 
stars 



INTRODUCTION 



The discovery of large numbers of extrasolar planets (exoplanets) has stimulated many 
observational and theoretical studies of their atmospheric chemical compositions. In particu- 
lar, the habitability of an exoplanet at a given distance from its host star is generally thought 
to depend on the total flux received from the host star, the availability of gre enhouse gasses, 
and the presence of biomarkers such as O2 and O3 (but see iFeng et al.l (120121 )). Estimates of 
the chemical composition of exoplanet atmospheres depend crucially on the near-ultraviolet 
(NUV, A =1700-3200 A), far-ultraviolet (FUV, A = 1170-1700 A), extreme ultraviolet (EUV, 
A =300-911 A), XUV (A =100-300 A), and X-ray (A < 100 A) radiation from the host star, 
which control important molecular photodissociation and photoionization processes. 

Early interest concerning the influence of a star on its substellar companions was focused 
on the erosion of volatile gasses from the primordial atmospheres of the terrestrial planets 
of our Solar System (especially Mars and Venus) by dissociative recombination powered by 
solar ultraviolet radiation and pickup-ion stripping by the solar wind. These processes would 
have been e normously enhanced wh en the Sun was young and more magnetically active than 
today (e.g. JZahnle fc Walkerlll982l ). Using a data base of solar- type sta rs in young cl usters of 
known age and extrapolating from solar fluxes over the sunspot cycle, lAyresl (119971 ) showed 
that photoionization rates relevant to the Martian situation scale as power laws in time 
(as anticipated in the earlier Zahnle & Walker paper). Furthermore, the UV flux behavior 
could be reliably traced back to the epoch of the Late Heavy Bombardment (at solar age 
~800 Myr), because by then rotation rates of you ng G stars had funnele d into a narrow 
distribution from an initially more diverse behavior ( IDelorme et al.ll2011(@l ). It is the stellar 
spin that largely governs the degree of magnetic activity, and the decay of that spin by wind 
braking that causes the rapid decline of the UV and X-ray activity with increasing age. 



More recently, the "Sun in Time" project (IRibas et al.ll2005l . I201CH ) and lSanz-Forcada et al 
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( 120 111 ) have extended the earlier work to a wider range of high energy measurements, es- 
pecially including the important 912-1170 A Far Ultraviolet Spectroscopic Explorer (FUSE) 
band, and to a broader range of spectral types (GO to G5), utilizing field stars whose ages 
must be inferred indirectly, but are more accessible to observation than the mostly distant 
cluster members. These studies also showed that in different wavelength bands, power laws 
could be used to describe the decay of chromospheric and coronal radiation over time. 



Rib as et al.l ( 20051 ) called attention to the importance of the H I Lyman-q line (1215.67 A), 



which is by far the brightest FUV emission line. They showed that for the Sun, the Lyman-a; 
line flux is about 20% of the total flux between 1 and 1700 A. The relative importance of 
the Lyman-a line is even more important for cooler stars as the photospheric emission at 
A > 1700 A decreases rapidly with decreasing effective temperature. 

Since important molecules in planetary atmospheres, including H 2 0, CO2, and CH 4 , 
have photodissociation cross sections that peak below 1700 A0, Lyman-a radiation can be 
the dominant cause of photodissociation for these and other molecules. For example, the 
solar Lyman-a line is responsible for 51% of the photodissociation rate of H 2 by solar UV 
radiation between 1148 and 1940 A, and 71% of the total photodissociation rate of CH 4 by 
solar radiation between 56 and 1520 A. However, recent m odels of exoplanet atmospheres 
that include photoionization, e.g., models for WASP-12b (IKopparapu. Kasting. fc Zahnle 
20121 ). typically do not include realistic values for the stellar Lyman-q flu x. In this recent 
model, the FUV fluxes were obtained from the compilation of lPickled (119981 ) that is based on 
International Ultraviolet Explorer (IUE) spectra without a large correction for interstellar 
absorption in the Lyman-a line or geocoronal Lyman-a emission. New ultraviolet spectra 
of six M dwarf host stars observed with spectrographs on the Hubble Space Telescope by 
France et al.l (j2012bl ) provide an observational basis for more realistic photochemical models 
of exoplanet atmospheres. 

While exoplanet atmospheres absorb the Lyman-a flux from their host stars without 
attenuation, neutral hydrogen in the interstellar medium scatters most of the Lyman-a 
flux out of the line-of-sight between the star and Earth. One must therefore correct for 
interstellar absorption by recon structing the s tellar Lym an-q profile and fl ux. For pre- 



main sequence stars with disks, iHerczeg et al.l ( 120041 ) and ISchindhelm et al.l (120121 ) showed 
that the fluorescent H 2 emission lines pumped by Lyman-q provide a useful diagnostic for 
reconstructing the intrin sic Lyman-q emiss ion line. For 62 main-sequence and giant stars 
without H 2 fluorescence, IWood et al.l ( 120051 ) reconstructed Lyman-q fluxes using interstellar 
H I column densities and velocities inferred from the deuterium Lyman-q and metal lines. 



iMPI-Mainz-UV-VIS Spectral Atlas of Gaseous Molecules, 



http://www.atmosphere.mpg.de/enid/2295 
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France et al.l (I2012al ) developed an alternative reconstruction technique in which the widths 
and strengths of one or two Gaussian emission lines representing Lyman-a and the velocities 
and column densities of interstellar absorption features are iteratively varied to obtain a 
solution that best fits the observed Lyman-a line wings. Both techniques require high- 
resolution spectra with good signal/noise and minimal contamination by geocorona Lyman-a 
emission. The only available spectrograph that can obtain such data is the Space Telescope 
Imaging Spectrograph (STIS) on HST. High instrumental sensitivity, especially important 
for faint M dwarfs, increasing interstellar hydrogen column densities with distance, and 
diflicult-to-obtain HST observing time limit the number of stars for which Lyman-a fluxes 
can be reconstructed by these two techniques. 

The objective of this study is to find different methods for infering the Lyman-a flux 
incident on the atmospheres of exoplanets for stars when high-resolution Lyman-a spectra 
are not available. These methods should be applicable to a broad range of stars including 
M dwarfs, which are very numerous and could support habitable planets with short orbital 
periods. One proposed method is to use correlations of reconstructed Lyman-a line fluxes 
with other stellar features (lines of O I and C II) formed in the same temperature range 
as Lyman-a, lines of Mg II and Ca II formed at slightly lower temperatures, and lines of 
C IV formed at somewhat higher temperatures. Our hypothesis is that ratios of the Lyman- 
a line flux to fluxes in these other lines should be similar for stars in limited ranges of 
spectral type with only a small dependence on stellar activity as indicated by the emission 
line fluxes. Empiri cal support f or this hypothesis comes from the similar power law slopes of 
UV em ission lines (Avresj 1997 ; Ribas et al. 2005 ) and the correlation of Lyman-a and Mg II 



fluxes ( Wood et al. 20051 ). The similar shapes of the chromosphere and tr ansition region 



therm al structures for regions on the Sun with very different heating rates ( jFontenla et al. 



20111 ) provides theoretical support for our hypothesis. 



The paper is structured as follows. In section 2 we describe our selection of targets 
that have reconstructed Lyman-a and other UV emission line fluxes. Section 3 describes 
our procedure for reconstructing the intrinsic Lyman-a flux using correlations with other 
UV emission lines and the errors and uncertainties of this method. Section 4 describes an 
extension of this procedure using correlations with ground-based Ca II H and K line fluxes, 
and section 5 describes how correlations with X-ray fluxes can be used to predict the intrinsic 
Lyman-a flux. Section 6 describes how the stellar effective temperature and rotation rate 
can provide rough estimates of the intrinsic Lyman-a flux when no emission line fluxes are 
available. The final section lists our conclusions. In a subsequent paper, we will apply 
the flux correlation technique to estimating stellar EUV radiation, which is responsible for 
photoionizing hydrogen and other species, heating outer atmospheres, and driving mass loss 
from exoplanets. 
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TARGET SELECTION 



Our target list consists of all F5 to M5 main se quence stars for which reconstructed 
Lyman-a fluxes are now available. IWood et al.l ( 120051 ) reconstructed Lyman-a line profiles 
observed by STIS in its medium (spectral resolution « A/45,000) and high-resolution (« 
A/100,000) echelle modes using interstellar deuterium Lyman-a and metal absorption lines 
to model the interstellar H I column density as a function of wavelength across the hydrogen 
Lyman-a line. They provided reconstructed Lyman-a line fluxes for 40 main-sequence stars 
with spectral types F5 to M5, together with Mg II h and k line and X-ray fluxes for most 
of these stars. We also include five M dwarf stars that are known to host exoplanets. The 
Lyman-a fluxes for these stars were reconstructed from STIS grating and echelle spectra 
using an it erative technique to measure the intrinsic Lyman-a profile and the interstellar 
absorption (IFrance et al.ll2012al ). 



For most these stars, we have measured fluxes of the H I 1215.67 A, O I 1304.86 A + 
1306.03 A, C II 1335.71 A, and C IV 1548.19 A + 1550.77 A lines. In the semi-emp irical solar 
chromosphere and transition region model computed by lAvrett fc Loeserl (120081 ). the peak 
contributions to the line center emission for these lines are at temperatures near 40,000 K, 
6,700 K, 29,500 K, and 68,000 K, respectively. However, the wings of these lines are formed 
over a range of cooler temperatures. We did not use fluxes of the O I 1302.17 A and C II 
1334.53 A lines as these lines show strong interstellar absorptio n. We extracted th e O I, C II 
and C IV line fluxes from Cosmic Origin Spectrograph (COS) (IGreen et al.ll2012l ) data with 
spectral resolution (« A/17, 000) available through the Mikulski Ar chive for Spa ce Telescopes 
(MAST)@ and the well-calibrated STIS spectra from the Star CAT (Uyreslboioj ) websitfl In 
a few cases we used data obtained with the Goddard High Resolution Spectrograph (GHRS) 
instrument on HST. The O I data obtained with COS are not usable because of airglow 
contamination. We subtracted the underlying co ntinuum from the em ission line fluxes. Line 
fluxes for the solar-mass stars were measured by iLinsky et al.l (120121 ). 



We also include solar spectral irradiance data observed with the Solar Radiation and 
Climate Experi ment (SORCE) on the Solar-Stellar Irradiance Comparison Experiment II 
(SOLSTICE II) JWoods et al.lbood : ISnow et al.ll2005h . These are integrated sunlight spectra 
directly comparable with the stellar data. The 2008 April 10-16 data refer to the very quiet 
Sun, and the 2003 October 27 and November 10 data refer to the Sun when it was very 
active. 



2 http://archive. stsci.edu 

3 casa. colorado.edu/~ayres/StarC AT/ 
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Table 1 summarizes the stellar and emission line fluxes f(line) in ergs cm -2 s" 1 at a 
distance of 1 AU for 45 stars, the quiet Sun, and the active Sun at two different times. Since 
most of the stars are variable, especially in the FUV, the fluxes refer to a single time and are 
not time averages. For most of the stars the STIS Lyman-a spectra and the COS spectra of 
the O I, C II, and C IV were obtained at different times. Since stellar UV radiation varies 
with time especially for M dwarfs and active warmer stars, ratios of the Lyman-a flux to 
other emission lines will have systematic errors compared to the stars for which all of the 
emission lines were observed at nearly the same time. 

3. CORRELATIONS OF RECONSTRUCTED LYMAN-a FLUX WITH 
EMISSION LINES FORMED IN THE CHROMOSPHERE AND 

TRANSITION REGION 

3.1. F5 V to G9 V stars 

Table 1 lists the measured fluxes for 21 F5-G9 stars and the Sun at three different 
times and activity levels. We plot in Figure 1 the flux ratio R(C IV) = f(Lyman-a)/f(C IV) 
vs. f(C IV) with line fluxes measured in ergs cm" 2 s" 1 at a distance of 1 AU. Figures 2-5 
show similar plots of R(line) = f(Lyman-a)/f(line) vs. f(line) using the previously described 
C II, O I, and Mg II lines. The R(line) ratios for the F5-G9 stars in Figures 1-4 follow 
tight trajectories of decreasing R(line) with increasing f(line). In Table 1 we list the iron 
abundances relative to hydrogen [Fe/H] listed in the SIMBADj data base. Three of the F5- 
G9 stars (HR 4657, \ Her, and r Cet) have large iron depletions (defined as [Fe/H] < 0.3), 
and often have higher R(line) values than other stars with similar fluxes and [Fe/H] values 
close to the solar. This behavior is expected since R(line) should increase as the line's metal 
abundance decreases, but whether or not R(line) is inversely proportional to [Fe/H] will be 
considered in the next section. We therefore exclude the low [Fe/H] stars when computing 
least-squares fits for the remaining F5-G9 stars. The coefficients for these fits, log[R(line)] 
= A + Blog[f(line)], are listed in Table 2. We note that the least-squares fits, which do not 
include the solar data, are close to the solar data. Also, the a Cen A ratios are similar to the 
quiet and active Sun ratios even though the solar and stellar data are measured by different 
instruments. Table 2 shows that the mean dispersions of the F5-G9 stars about the best fit 
lines are 18-24%. 



4 lit t p : / / simbad . u- st r asbg . fr / simbad 
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3.2. KO V to K5 V stars 

There are a total of 16 stars in this group: 15 have Mg II fluxes, and 8 have C IV, C II 
fluxes and O I fluxes. The least-squares fits to the R(line) ratios for the KO V - K5 V stars 
have similar slopes to the fits to the R(line) ratios for the F5-G9 stars and similar dispersions 
(15-30%) about the fit lines. The KO V-K5 V stars have a wide range in rotational periods 
(0.38 days for Speedy Mic to 43 days for a Cen B), but there is no strong trend of R(line) 
with rotational period and thus with activity. Eight of these stars are classified in SIMBAD 
as either BY Dra type varables or RS CVn type variables, but we can find no information 
on the companion stars which are likely fainter and cooler than the primary stars. 



3.3. MO V to M5 V stars 

We analyze spectra of nine MO V to M5.5 V stars with reconstructed Lyman-a fluxes. 
Five of these are exoplanet host stars. We plot R(line) vs. log[f(line)] in Figures 1-3 and 
5. The dispersion of the data about the fit lines is much larger for the M stars than for the 
F5-K5 stars primarily due to their much larger time variability as described in the following 
section. The dispersion in R(0 I) for the M stars, however, may be unrepresentative given 
the small number of data points. 



3.4. Errors and Uncertainties 



In this analysis, we identify four causes of uncertainty in the R(line) values: flux mea- 
surement errors, Lyman-a line reconstruction errors, uncertain atomic abundances, and un- 
known time variability. We can estimate the error magnitudes or at least identify the first 
three types of errors, but time variability errors cannot be measured with the existing data 
and are therefore systematic. 

The smallest errors are typically flux-measurement errors for the C IV, C II, O I, and 
Mg II lines. These are typically less than 2% for the STIS data in StarCAT and the COS 
data sets, but can be as large as 5-20% for the faintest M dwarfs. Uncertainties in the re- 
construction of the Lyman-a lines are usually larger than the line flux measurement errors. 
Wood et al.l (120051 ) estimated 20% for typical errors in reconstructing Lyman-a line fluxes. 
In order to compare the tw o reconstruct i on m ethods, we re-fit the observed Lyman-a line 
profile s of several stars in t he lWood et al.l (120051 ) sample with the iterative least-squares tech- 
nique ([France et al.ll2012al ). For the favorable case of AU Mic, a star with high signal/noise 
data and a single ISM absorption feature, the agreement in the reconstructed fluxes between 
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the two techniques is within 5%. For the unfavorable case of AD Leo, a star with several ISM 
velocity components along the line of sight, the disagreement is about 3 0%. In Figures 1-5 , 
we plot representative errors of 20% for the R(line) ratios based on the IWood et al.l (120051 ) 
reconstruction, or 30% for the faint M dwarfs based on the lFrance et al.l (J2012a]) reconstruc- 
tion. As shown Table 2, the mean dispersions of the F5-G9 and K0-K5 stars (ignoring the 
low metal abundance stars) are in the range 18-30%, which is consistent with our estimates 
of the Lyman-a reconstruction errors and smaller time variability errors. 

We have excluded stars with low iron abundances ([Fe/H]< 0.3) from the sample of 
stars used to make the least-squares fits. The [Fe/H] abundances for the F, G, and K 
stars were obtained from the SIMBAD data base. Accurate abundance measurements for 
M dwarfs require special techniques given the complexity of their spectra. Recent work in- 
cludes calibration of th e offset of stars in the H-R diagram relative to the main sequence 
( Johnson fe Apps 2009|) and m easurements of Na I and C a I absorption lines in near-IR spec- 



tra (IRojas-Ayala et al.ll2010l ). iRojas-Ayala et al.l (120121 ) have revised their earlier work and 
compared their new results to earlier studies. We therefore adopt their [Fe/H] measurements 
for AD Leo, EV Lac, and the exoplanet h ost stars GJ 876, GJ 58 1, and GJ 436. We adopt 



the [F e/H] value for GJ 832 measured by ( IRojas-Ayala et al.ll2010l ). ICayrel de Strobel et al. 



( 1200 ll ) lists [Fe/H] =-0.54 for GJ 667C, but we do not expect this star to have a significant 
under abundance as the other host stars show normal or superrich abundances. The young 
star AU Mic should have near solar abundances, and Proxima Cen should have the near-solar 
abundances of a Cen A and B. 

For each low iron abundance star, Table 4 lists the value of [Fe/H], and for each spectral 
line the log of the difference, which we refer to as "delta", between R(line) and the least- 
squares fit line at the value of f(line). If the deltas were only due to low metal abundance 
and the abundances of carbon, oxygen, and magnesium relative to the Sun were the same as 
for iron, then the observed deltas should be the negative of the log iron depletions. This is 
not the case, but the data in Table 4 show some interesting trends. For example, the deltas 
seen in the C IV and C II lines (and to some extent in the O I line) are nearly the same for 
each star, implying that the abundances of carbon and oxygen may be important factors in 
determining the deltas. However, the deltas are much smaller than -[Fe/H] for the F5-G9 
and K0-K5 stars. Thus metal depletion cannot be the only factor determining the deltas, 
and the thermal structures of the atmospheres of stars of different spectral type likely play a 
role. Since the many lines of Fe, Mg, an d Ca are the main cooling agents of the lower solar 
chromosphere (lAnderson fe Athaylll989l ). lower abundances of these elements will change the 
energy balance and thus the thermal structure. 



Since the deltas for the Mg II lines of the F5-K5 stars are mostly smaller than those 
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for the C IV, C II, and I lines, there will likely be a smaller uncertainty when estimating 
R(Mg II) from the least-squares fits compared to using R(line) for the C IV, C II, and O I 
lines. However, it is important to correct for interstellar absorption in the Mg II lines, which 
requires high-resolution spectra. Since interstellar absorption is minimal for the C IV lines, 
the C II 1335 A line, and the O I 1304 and 1306 A lines, lower resolution spectra of these 
lines can be used without interstellar corrections for estimating R(line) and thus the intrinsic 
Lyman-a flux. 

In most cases, the STIS spectra containing the Lyman-a line and the COS or STIS 
spectra used for extracting the other line fluxes were obtained at different times. For these 
stars the R(line) ratios will depend on time variations in the stellar UV emission. The 
importance of time variability for different types of stars can be estimated from the mean 
dispersion and RMS dispersions of the R(line) ratios about the least-squares fits. These 
dispersions for the M stars listed in Table 2 are much larger than for the F5-K5 stars and 
much larger than the 20-30% uncertainty associated with the Lyman-a reconstructions. This 
behavior is expected since M stars are more variable and flare more often than the warmer 
stars. 

4. CORRELATION OF RECONSTRUCTED LYMAN-a FLUX WITH Ca II 

EMISSION 

The hydrogen Ha (6563 A) and Ca II H and K (3933 and 3968 A) lines formed in the 
chromosphere are observable by ground-based telescopes. However, the Ha line is difficult 
to analyze since increasing heating first deepens the absorption line and then produces an 
emission feature. We consider instead the Ca II lines because the emission in the centers of 
the broad photospheric absorption lines is a good indicator of the chromospheric heating rate. 
An important difficulty is distinguishing the chromospheric emission inside of the HI and Kl 
features that define the emission core from the photospheric emission that would be present 
in the absence of a chromosphere. For active stars and especially M dwarfs, this is not a 
large uncertainty as the photospheric emission is faint compared to the bright chromospheric 
emission. For the less active G stars, the photospheric emission is comparable to or larger 
than the chromospheric emission. It is therefore essential to correct for the photospheric 
emission as we wish to correlate chromospheric Ca II emission with the intrinsic Lyman-a 
flux. 



Different authors has estimated the photospheric emission in and near t he co re of 
the Ca II lines in di fferent ways. Using high-resolution spectra, iLinsky et al.l (Il979l ) and 
Pasquini et al.l (119881 ) measured the flux in the Ca II emission line cores between the HI 
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and Kl minima features and subtracted the small amount of flux in these wavelengt h in- 
tervals predicted by r adiative equilibrium model photospheres. iRobinson et al.l (Il990h and 
Browning et al.l ( 120101 ) fitted the photospheric H and K absorption lines with Gaussians and 
then subtracted the flux interpolated in the line cores from the measured Ca II emission. 
This approach somewhat overcorrects for the photospheric emission, but the error is small 
for the very cool stars that they observed. The resulting chromospheric Ca II fluxes at 1 AU 
are listed in columns 6 and 7 of Table 3. 

Other observers have measured the S index, which is the flux in 1.09 A bandpasses 
centered on the H and K line emission cores divided by the flux in continuum windows. 
Hartmann et al.l (119841 ) provided a correction for the photospheric flux within the bandpass 
but outside of the HI and Kl features using high-resolution spectra of stars with weak H 
and K emission. Colu mns 8 and 9 of Ta ble 3 show the time-avera ged Ca II fluxes ob- 
tained at Mt. Wilson bv lNoves et al.l (119841) a nd iBaliunas et al.l (119951 ) after subtracting the 
photospheri c flux following lHartmann et al.l (I1984J) . Co lumns 10 and 11 show Ca II fluxes 



Ca II), using the avera ge of the Ca II fluxes for all of the data in Table 3, except 



Cincunegui et al.l ( 120071 ) fluxes which appear to be systematically smaller than the 



obtained by iHenrv et al.l (119961 ). lHall et al.l (120071 |2009|), and lLockwood et al.l (12007D using 
a similar ap proach. The last c o lumn in Table 3 lists the fluxes of Southern hemisphere stars 
obtained by ICincunegui et al.l (120071 ) . We plot in Figure 6 the ratio of Lyman-a to Ca II 
flux, R 
for the 

fluxes obtained from the other sources.. 

The R(Ca II) data in Figure 6 for each spectral type bin show relatively small scatter 
about the least-squares fit lines. We find that the dispersion (see Table 2) is smallest for the 
K stars (only 20.6%), and larger for the variable M stars. The dispersion of 42% for the G 
stars likely results from uncertainty in the photospheric flux correction which is large for G 
stars. Thus the Ca II H and K line flux can be used to estimate the intrinsic Lyman-a flux 
from the fit lines, provided one can remove the photospheric component from the observed 
line core fluxes. 



5. CORRELATION OF RECONSTRUCTED LYMAN-a FLUX WITH 

X-RAY EMISSION 

Except for the highly variable M dwarfs for which the Lyman-a and other emission 
lines were observed at different times, the correlations of reconstructed Lyman-a flux with 
emission lines formed at similar or somewhat higher temperatures in the chromosphere and 
transition region generally have small scatter about the least-squares fit lines. This behavior 
results from increasing mechanical heating shifting the chromospheric thermal structure 
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to higher densities while keeping a similar shape (IFontenla et al.l 1201 ll ). Thus al l of th e 
emission lines brig hten together roughly proportional to density squared. lAyresI ( 119971 ). 



Ribas et al.l ( 120051 ). and others have noted that slopes of power- law relations between X- 
ray and chromosphere line flux are generally near 2.0 rather than a linear relation. Also, 
soft X-ray emission as measured, for example, by ROSAT is highly sensitive to coronal 
temperature and is more time variable than chromospheric emission. For these reasons, we 
did not anti c ipate a tight correlation between reconstructed Lyman-a flux and X-ray flux. 
Wood et al.l (120051 ) plotted X-ray flux vs reconstructed Lyman-a flux for F-M dwarfs and 
giants showing correlations, but with large scatter for each type of star. 

We plot in Figure 7 the R(X)=f(Lyman-a)/f (X-ray) ratios vs. X-ray flux at 1 AU for 
the stars in Table 1 using the ROSAT PSPC X-ray fluxes cited by IWood et al.l (120051 ) or 
more recent X-ray flux measurements using Chandra and XMM-Newton (see the references 
cited in the footnotes of Table 1). The X-ray flux for GJ 667C, which was obtained with the 
HRI instrument on ROSAT rather than the PSPC, is uncerta in due to the close proximetry 
of GJ 667A and GJ 667B. lJudge. Solomon, fc Ayred (120031 ) estimated solar X-ray fluxes 
through the ROSAT PSPC bandpass from full disk solar observations by the SNOE-SXP 
instrument. For the quiet Sun, moderately active Sun, and active Sun, we use their X-ray 
luminosities, logL^ = 26.8, 27.75, and 27.9, respectively. 

The Lyman-a/X-ray flux data show correlations with small scatter for the F5-G9 and 
K0-K5 dwarfs. The quiet and active solar data are consistent with the least-squares fit for 
the other F5-G9 stars. The mean dispersions listed in Table 2 are 29.8% for the F5-G9 
dwarfs and 22.2% for the K0-K5 dwarfs, which are similar to the mean dispersions for the 
C IV, C II, O I, Mg II, and Ca II lines. The mean dispersion for the M dwarfs is much larger 
than for the F5-K5 stars, although the slope of the fit line is similar to that of the warmer 
stars. X-ray time variability is the most likely cause of the large scatter of the M stars as 
the X-ray and Lyman-a data were obtained at different times. Despite the factor of 2-3 
scatter about the M star fit line, the similar slope of the fit line to those for the warmer stars 
suggests that one can use the M dwarf fit line to estimate the Lyman-a flux of an M dwarf 
at the time of an X-ray measurement with an uncetainty much smaller than the dispersion 
shown in Table 2. We note that GJ832 is plotted very close to the quiet Sun in Figure 7, and 
that GJ832 (Ml. 5 V) and AU Mic (M0 V) are consistent with the fit lines for the K0-K5 
stars. This suggests that the R(X) vs. X-ray flux correlation for K stars may also be useful 
for estimating the intrinsic Lyman-a flux of early M dwarfs. 
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6. CORRELATION OF THE RECONSTRUCTED LYMAN-a FLUX WITH 
STELLAR EFFECTIVE TEMPERATURE AND ROTATION 

We now consider whether the stellar Lyman-a flux of main- sequence stars can be esti- 
mated simply from the stellar effective temperature (T e g). We plot in Figure 8 the Lyman-a 
flux for all of the stars in Table 1 vs. T e s. At each value of T e ff, there is a dispersion in the 
Lyman-a flux that increases from a factor of 4 near T e g = 6000 K to a factor of 1000 for 
the M stars. The range in Lyman-a flux at a given T e g is reduced significantly by grouping 
the stars according to their rotation period (P r ot)j which is a rough measure of magnetic 
heating rates in the chromosphere and corona. We show least squares fits, log f(Lyman-a) 
= A + BT ef f, for the stars in three groups: fast rotators (P rot = 3-10 days), moderate ro- 
tators (10-25 days), and slow rotators (> 25 days). The coefficients for these fits are listed 
in Table 5. Compared to the least-squares fits, the Lyman-a flux dispersion within each 
rotational period group is 32-85% and largest for the M stars. Therefore, estimates of the 
Lyman-a flux based only on a star's T eff and P rot are reasonably accurate for G stars but 
not for the cooler stars. 

In Figure 9, we plot the Lyman-a flux in the habitable zone (HZ) of an exoplanet, where 
the distance from the star to the habitable zone is estimated approximently as proportional 
to the stellar luminosity, d H z = (-R*/-R©)(^*/^©) 2 AU. The location of the HZ, the region 
around a star where life forms are thought to be possible, involves many other considerations 
including greenhouse gases, orbital eccentricity and stability, and atmospheric chemical com- 



position (IKasting Sz Catling! 120031 ). In this HZ plot, the dispersion of the Lyman-a flux is 
reduced somewhat for the M stars (see Table 5). Lyman-a fluxes in the HZ trend significantly 
higher with decreasing T e g, and the quiet Sun has nearly the smallest value of Lyman-a flux. 
Exoplanets of most of the stars discussed in this paper receive significantly larger Lyman-a 
fluxes and thus have higher molecular dissociation rates in their atmospheres. This is espe- 
cially true for M dwarfs, which have smaller radii and lower effective temperatures than the 
Sun. Exoplanets in the HZs of the five M dwarfs in this study receive about 10 times the 
Lyman-a flux as the Earth receives from the quiet Sun. We therefore conclude that large 
FUV fluxes are received by exoplanets in their HZs from their M dwarf host stars contrary 
to some previous assumptions. Also, the Lyman-a fluxes for stars with exoplanets do not 
appear to differ from stars without discovered exoplanets in these two plots. 



7. CONCLUSIONS 



We have identified five methods for reconstructing or estimating the intrinsic Lyman-a 
stellar fluxes of main-sequence stars between spectral types F5 and M5. Whether a given 
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method can be used and its accuracy de pends on the qualit y of the available data. The first 



and most accurate method developed by lWood et al.l (120051 ) requires high-resolution spectra 



of the Lyman-a line and knowledge of the interstellar veloc ity structure based on high- 



resolution spectra of the deuterium Lyman-a and metal lines. I Wood et al.l (120051 ) estimates 
typical errors of 20% in the reconstructed Lyman-a line fluxes using this method, but at 
present only the STIS instrument on HST can provide such data for nearby stars. 



The second method developed by France et al.l (j2012al ) also requires high-resolution 



spectra of the hydrogen Lyman-a line, but does not require spectra of the deuterium Lyman- 
a line or any other interstellar absorption line. This technique solves for both the intrinsic 
Lyman-a line parameters and the interstellar absorption simulataneously. When the inter- 
stellar absorption has only one velocity component, the technique can be as accurate as the 
first, but it fails when the interstellar velocity structure has many components, which is not 
known in the absence of high-resolution interstellar absorption lines. 

In this paper, we considered three additional methods that can be used with different 
accuracy when there is no available high-resolution spectrum of the Lyman-a line to serve 
as the basis for reconstruction. The third method described in Sections 3 and 4 requires flux 
measurements of the stellar C IV, C II, O I, Mg II, or Ca II lines and the best fit correlations 
of these lines with fluxes of the reconstructed Lyman-a lines in our data set. This method 
estimates Lyman-a fluxes with 18-25% uncertainty for F5-K5 dwarf stars, provided one 
corrects high-resolution spectra of the Mg II lines for interstellar absorption. This method 
is based on the hypothesis that the ratios of the Lyman-a line flux to C IV, C II, O I, Mg II, 
and Ca II line fluxes, R(line), for stars of similar spectral type depend only gradually on line 
flux. We find that this hypothesis is valid for late-type dwarf stars with approximently solar 
abundances. Using this method for F5-K5 V stars, we find that the dispersions of R(C IV), 
R(C II), R(0 I), and R(Mg II) about the least-squares fits are consistent with the likely 
errors in reconstructing the intrinsic Lyman-a fluxes. The dispersion of R(Ca II) about the 
fit line for the F5-G9 stars is larger than for the other lines because of uncertain estimates 
of the photospheric flux. 

Most dispersions for the M0-M5 V stars are significantly larger, probably due to stellar 
variability between the time the Lyman-a line was observed with the STIS instrument and 
the other lines were observed with COS or STIS with a different grating setting. Even with 
this time variability, the fit lines for the O I, C II, and C IV lines should provide estimates of 
the intrinsic Lyman-a flux for a given star with an uncertainty less than a factor of two. We 
suggest that estimating the intrinsic Lyman-a flux from the Mg II line flux and the least- 
squares fit will lead to smaller uncertainty provided one has high-resolution Mg II spectra 
with which to estimate interstellar absorption in these lines. 



- 14 - 



It is important, however, to recognize the limitations of this correlation method. The 
errors associated with measuring fluxes of the Mg II, O I, C II, and C IV lines are generally 
small at the time of these line flux measurements, although geocoronal emission through 
the large COS aperture can be important for the O I lines. Since the method does not use 
direct measurements of the Lyman-a line, there is no reconstruction error. In addition, there 
is no time variability error as the Lyman-a flux is the scaled value at the time the other 
emission lines are observed. The largest uncertainty is the stellar metal abundance. Stars 
known to have low metal abundances generally have R(line) ratios above the scaling relation 
predictions. Since the effect appears to be smallest for the Mg II lines, correlations with the 
Mg II line fluxes may provide more accurate intrinsic Lyman-a line fluxes. Observations of 
more stars, especially M stars, are need to better understand the errors in R(line) associated 
with low metal abundances. 

When no UV or Ca II line fluxes are available but there are X-ray measurements with 
an energy range similar to that of the ROSAT PSPC, a fourth method can provide estimates 
of the intrinsic Lyman-a flux from least-squares fits to the intrinsic Lyman-a/X-ray flux 
ratio vs. X-ray flux. The mean dispersions about the fit lines for F5-G9 dwarfs and K0-K5 
dwarfs are 20-30%, but the mean dispersion for the M dwarfs is much larger again due to 
the large time variability of X-ray emission and the comparison of X-ray and Lyman-a data 
obtained at different times. Since the slope of the Lyman-a/X-ray flux ratio vs. X-ray flux 
for the M stars is similar to that obtained for the warmer stars, we suggest that one can 
obtain a factor of two estimate of the intrinsic Lyman-a line flux at the time of the X-ray 
measurement from the correlation fit line. 

Even when no Lyman-a, UV emission lines, or X-ray data are available, one can use a 
fifth method to estimate the intrinsic Lyman-a flux for F5-M5 stars based only on the star's 
effective temperature and some measure of stellar activity. The plot of intrinsic Lyman-a flux 
vs. stellar effective temperature shows more than an order of magnitude range in Lyman-a 
fluxes at a given effective temperature, but comparing data for stars of similar rotational 
period, a good measure of stellar activity, significant reduces the range. Comparison of 
Lyman-a fluxes vs. effective temperature for stars of similar rotational period as viewed 
from their habitable zones further reduces the mean dispersion about the fit lines to 30- 
40%. Thus even in the absence of any spectroscopic data, this fifth method can provide 
useful estimates of the intrinsic Lyman-a flux of F5-M5 dwarf stars. 

This work is supported by NASA through grants NNX08AC146, NAS5-98043, and HST- 
GO-11687.01-A to the University of Colorado at Boulder. We thank Martin Snow for pro- 
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Table 1. Stellar Line Fluxes (ergs cm 2 s 1 ) at 1 Astronomical Unit 



Star a 


HD 


[Fc/H] 


Spec Type b 


d(pc) b 


Lyman-a c 


Mg II h+k d 


O I 130.4+130.6 


C II 133.6 


C IV 154.8+150.1 


X-ray 


Ref ° 


Procyon 


61421 


-0.02 


F5 IV-V 


3.50 


77.1 


267 


1.77 


2.58 


4.62 


11.5 


5 


HR 4657 


106516 


-0.70 


F5 V/L 


22.6 


27.8 


62.6 


0.418: 


0.320: 


0.493: 


5.43 


1 


C Dor 


33262 


-0.15 


F7 V 


11.7 


46.5 


108 


0.627 


1.050 


1.958 


16.7 


2 


X Her 


142373 


-0.50 


F8 V/L 


15.9 


22.0 


45.2 


0.235 


0.115 


0.171 


0.312 


1 


— 


28033 


+0.11 


F8 V 


46.4 


24.7 


50.4 


— 


— 


— 


19.2 


— 


X 1 Ori 


39587 


-0.09 


GO V 


8.66 


41.6 


79.8 


0.474 


0.684 


1.247 


37.3 


2 


HR 4345 


97334 


-0.01 


GO V 


21.93 


42.8 


108 


0.569 


1.016 


1.631 


39.9 


2 


SAO 136111 


73350 


+0.07 


GO V 


23.98 


32.8 


65.8 


0.296 


0.429 


0.737 


19.3 


2 


V993 Tau 


28205 


+0.12 


GO V 


47.01 


55.5 


140 


0.901 


1.286 


2.46: 


41.4 


1 


V376 Peg* 


209458 


-0.06 


GO V 


49.63 


15.7 


— 


— 


— 


— 


— 


— 


Quiet Sun (4/2008) 




+0.00 


G2 V 


— 


5.95 


18.2 


0.0793 


0.0862 


0.129 


0.224 


3 


Active Sun (11/2003) 




+0.00 


G2 V 


~ 


7.04 


" 


0.0881 


0.1041 


0.152 


1.99 


3 


Active Sun (10/2003) 




+0.00 


G2 V 




9.15 




0.1107 


0.1497 


0.212 


2.85 


3 


a Cen A 


128620 


+0.25 


G2 V 


1.325 


7.54 


29.7 


0.1032 


0.1057 


0.161 


0.117 


4,13 


HR 2882 


59967 


-0.19 


G4 V 


21.82 


55.9 


90.2 


0.550 


0.821 


1.881 


41.9 


2 


61 Vir* 


115617 


+0.00 


G5 V 


8.555 


5.26 


14.1 


0.0488 


0.0427 


0.0603 


0.265 


i — 1 

loo 


K Cct 


20630 


+0.09 


G5 V 


9.14 


30.0 


53.0 


0.366 


0.556 


0.873 


25.6 


2 


HR 2225 


43162 


+0.00 


G5 V 


16.72 


41.0 




0.396 


0.479 


0.835 


48.1 


1 


HR 6748 


165185 


-0.06 


G5 V 


17.55 


48.9 


101.1 


0.495 


0.687 


1.123 


53.5 


2 


SAO 254993 


203244 


-0.21 


G5 V 


20.42 


43.8 


59.1 


0.311 


0.539 


0.940 


20.2 


1 


SAO 158720 


128987 


+0.01 


G6 V 


23.68 


34.4 


60.1 


0.296 


0.381 


0.510 


14.3 


1 


t Cet 


10700 


-0.43 


G8 V/L 


3.65 


5.66 


7.93 


0.0417 


0.0229 


0.0327 


0.176 


5 


£ Boo A 


131156A 


-0.13 


G8 V 


6.70 


35.3 


61.9 


0.274 


0.465 


0.815 


28.3 


1 


SAO 28753 


116956 


+0.03 


G9 IV-V 


21.9 


33.0 


57.5 


0.336 


0.504 


0.781 


24.7 


1 


a Cen B* 


128621 


+0.24 


K0 V 


1.255 


10.1 


19.1 


0.0809 


0.0925 


0.132 


0.533 


5,13 


DX Leo 


82443 


-0.23 


K0 V 


17.8 


31.1 


79.1 








59.7 




70 Oph A 


105341 


-0.08 


K0 V 


5.09 


23.6 


30.6 


0.222 


0.296 


0.431 


6.62 


5 


HR 8 


166 


+0.15 


K0 V 


13.67 


37.9 


57.0 


0.374 


0.592 


1.036 


33.0 


1 


e Eri* 


22049 


-0.08 


Kl V 


3.216 


21.5 


27.2 


0.145 


0.179 


0.274 


5.63 


1,9 


40 Eri A 


26965 


-0.27 


Kl V 


4.985 


7.33 


12.3 








1.15 




36 Oph A 


155886 


-0.39 


Kl V/L 


5.464 


18.0 


14.0 


0.0816 


0.0873 




3.72 




HR 1925 


37394 


+0.14 


Kl V 


12.28 


29.3 


45.0 


0.228 


0.277 


0.462 


14.4 


1 


_* 


189733 




Kl V 


19.25 


11.8 




0.202 


0.274 




5.34 


5,9 


EP Eri 


17925 


+0.08 


K2 V 


10.35 


27.6 


61.1 






0.917 


32.9 


5 


LQ Hya 


82558 




K2 V 


18.62 


59.1 


75.2 






4.744 


243. 


5 



Table 1 — Continued 





1 1 1 v 


f IA, / LI 1 


opec lype 




Lyman-ct 


Mn- TT Vi-Ll^d 
ivlg 11 ll^K 


f| T 1 Qf) A J_1 QD ft 


P IT 1 QQ (\ 

11 loo.u 


O IV IO^.OtIOU.1 


X-ray 


ixei 


V368 Cep 


220140 




K2 V 


19.20 


46.9 


86.6 








275. 




PW And 


1405 




K2 V 


21.9 


47.1 


51.8 






3.154 


187. 


5 


Speedy Mic 


197890 


-1.49 


K2 V/L 


52.2 


214 


190 


- 


- 




4255. 




61 Cyg A 


201091 


-0.35 


K5 V/L 


3.487 


8.90 


7.35 


0.0353 


0.0323 




0.498 


5,13 


e Ind 


209100 


-0.20 


K5 V 


3.622 


17.3 


8.43 








0.871 




AU Mic 


197481 




MO V 


9.91 


43.0 


17.6 


0.360 


0.525 


1.035 


70.9 


1,8 


AD Leo 


GJ 388 


+0.28 


M3.5 V 


4.695 


9.33 


2.19 


0.0644 


0.137 


0.376 


19.1 


1,10 


EV Lac 


GJ 873A 


-0.01 


M3.5 V 


5.122 


3.07 




0.0163 


0.0402 


0.1094 


19.5 


1,10 


Proxima Cen 


GJ 551C 




M5.5 V 


1.296 


0.301 




0.000408 


0.00162 


0.00719 


0.142 


1,11 


GJ 832* 


204961 


-0.12 


Ml. 5 V 


4.954 


5.17 


0.311 




0.00278 


0.00762 


0.214 


6,9 


GJ 667C* 






Ml. 5 V 


6.9 


1.54 


0.113 








0.263 


12 


GJ 876* 




+0.18 


M5.0 V 


4.689 


0.409 


0.0315 




0.00874 


0.0186 


0.0574 


7.9 


GJ 581* 




-0.10 


M2.5 V 


6.3 


0.513 


0.0360 




0.000820 


0.00304 




6 


GJ 436* 




+0.04 


M3 V 


10.3 


1.571 


0.1038 




0.00182 


0.00623 


0.0334 


6.9 



a Exoplanet host stars listed in the exoplanets.org data base are marked with a * symbol. 
b Data from SIMBAD with M star metal abundances from the sources listed in Section 3.4. 

c Reconstructed intrinsic Lyman-o line fluxes for most stars JWood et al.ll20051) and for the GJ stars llFrance et al.ll2012bl) . 

d Mg II h and k core emission corrected for interstellar absorption for most stars llWood et al.|[2005h and for GJ stars (this paper). 

^Sources of the O I, C II, and C IV data: (1) StarCAT llAyres! l2010h , (2) iLinskv et al] ||2012T|. (3) Martin Snow (private commu nication) . (4) IPaeano et ah] 
l l2004h , (5) This p aper, (6llFrance et ahl [|2012al) , (71 iFrance et al.l I 2012'bVf8') X-ray fl ux in lSchneider fc Schmitj j20ld) . (9) X-ray flux in lSanz-Forcada et al.l ^Olll l. 
(10) X-ray flux in Iftobrade &; Schmitj 1120051 1. (11) X-ray flux ir 
iRobrade. Schmitt fc Favatal ||20L2TI 



Giidel et al.l ll2004T ). (12) X-ray flux from J. Schmitt (private communication) (13) X-ray flux in 



Table 2. Least-squares Fits to Line Flux Ratios 
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Flux a 


A b 


B b 


Mean 


RMS 


Group 
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Dispersion (%) 
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O I 


0.049 to 1.77 


1.872 


-0.198 


22.5 


27.2 


KO V - K5 V 


6 


O I 


0.08 to 0.37 


1.886 


-0.198 


22.3 


32.8 


MO V - M5 V 


4 


O I 


0.0004 to 0.36 


1.871 


-0.278 


16.9 


17.8 


F5 V - G9 V 


16 


Mg II 


14.1 to 267 


-0.291 


-0.0208 


24.2 


31.2 


KO V - K5 V 


12 


Mg II 


8.4 to 86.6 


0.338 


-0.318 


29.8 


39.3 


MO V - M5 V 


7 


Mg II 


0.032 to 17.6 


0.814 


-0.296 


27.8 


33.8 


F5 V - G9 V 


13 


Ca II 


11.4 to 149. 


-0.223 


0.080 


42.2 


49.6 


KO V - K5 V 


10 


Ca II 


5.9 to 39.4 


0.605 


-0.433 


20.6 


30.7 


MO V - M5 V 


7 


Ca II 


0.0076 to 6.53 


1.028 


-0.312 


39.9 


49.6 


F5 V - G9 V 


20 


X-ray 


0.12 to 52.5 


1.156 


-0.684 


29.8 


42.3 


KO V - K5 V 


15 


X-ray 


0.50 to 3090 


1.058 


-0.707 


22.2 


27.4 


MO V - M5 V 


8 


X-ray 


0.033 to 70.9 


0.431 


-0.573 


122. 


155. 



a Flux range of the C IV, C II, O I, and Mg II lines and X-ray flux in ergs cm 2 s 1 at a distance of 1 AU. 
b east-squares fit to log[f(Lyman-o)/f(line)] = A + B log[f(line)], where the line is C IV, C II, O I, or Mg II. 



Table 3. Ca II H and K Line Fluxes (ergs cm 2 s 1 ) at 1 Astronomical Unit 



Star a 


HD 


[Fe/H] 


Spec Type b 


Lyman-a c 


Ca II d 


Ca IP 


Ca II f 


Ca II s 


Ca II h 


Ca IP 


Ca IP 


Procyon 


61421 


-0.02 


F5 IV-V 


77.1 








149. 








HR 4657 


106516 


-0.70 


F5 V/L 


27.8 






46.4 


40.1 








X Her 


142373 


-0.50 


F8 V/L 


22.0 






31.5 


34.9 




25.5 




X On 


39587 


-0.09 


GO V 


41.6 






48.3 


53.0 




52.7 




HR 4345 


97334 


-0.01 


GO V 


42.8 






49.6 


54.5 




52.0 




Mean Sun 




+0.00 


G2 V 


6.5 


10.8 






16.1 




17.2 




a Cen A 


128620 


+0.25 


G2 V 


7.54 


16.6 










15.8 


7.75 


HR 2882 


59967 


-0.19 


G4 V 


55.9 










44.7 




27.4 


61 Vir* 


115617 


+0.00 


G5 V 


5.26 








11.3 


11.6 






K Cet 


20630 


+0.09 


G5 V 


30.0 


64.1 




39.4 


45.2 




43.0 




HR 6748 


165185 


-0.06 


G5 V 


48.9 










45.6 




21.8 


SAO 254993 


203244 


-0.21 


G5 V 


43.8 










33.0 




23.8 


r Cet 


10700 


—0.43 


G8 V/L 


5.66 






6.61 


7.24 


6.52 


6.46 


4.61 


£ Boo A 


131156A 


-0.13 


G8 V 


35.3 


31.5 




23.9 


26.9 


27.1 


23.3 




a Cen B* 


128621 


+0.24 


K0 V 


10.1 


5.75 








6.02 




3.05 


DX Leo 


82443 


-0.23 


K0 V 


31.1 








39.4 








70 Oph A 


105341 


-0.08 


K0 V 


23.6 


23.1 














e Eri* 


22049 


-0.08 


Kl V 


21.5 


19.7 




15.9 


16.8 


14.8 


15.9 


9.37 


40 Eri A 


26965 


-0.27 


Kl V 


7.33 


7.99 




7.42 


8.08 


6.41 




3.70 


36 Oph A 


155886 


-0.39 


Kl V/L 


18.0 






10.4 


10.5 


9.39 






HR 1925 


37394 


+0.14 


Kl V 


29.3 








23.4 








EP Eri 


17925 


+0.08 


K2 V 


27.6 






34.2 


34.1 


32.9 




22.2 


61 Cyg A 


201091 


-0.35 


K5 V/L 


8.90 


1.61 




3.0 


3.28 


7.89 


3.29 




e Ind 


209100 


-0.20 


K5 V 


17.3 


3.29 








8.13 




2.92 


AU Mic 


197481 




M0 V 


43.0 


8.60 


4.46 












AD Leo 


GJ 388 


+0.28 


M3.5 V 


9.33 




0.264 










0.958 


EV Lac 


GJ 873A 


-0.01 


M3.5 V 


3.07 




0.314 












Proxima Cen 


GJ 551C 




M5.5 V 


0.301 














0.0067 


GJ 667C* 






Ml. 5 V 


1.54 




0.103 












GJ 876* 




+0.18 


M5.0 V 


0.409 




0.00765 












GJ 581* 




-0.10 


M2.5 V 


0.513 




0.00828 












GJ 436* 




+0.04 


M3 V 


1.571 




0.0975 













a Exoplanet host stars listed in the exoplanets.org data base are marked with a * symbol. 
b Data from SIMBAD with M star metal abundances from the sources listed in Section 3.4. 

c Reconstructed intrinsic Lyman-ce line fluxes for most stars jWood et alj|2005h and for the GJ stars jPrance et alj|2012bh . 
d Ca II H and K chromospheric flux at 1 AU llLinskv et al.|[l979h ■ IPasquini et al. | {l98Sft and llRobinson et al.lll990h . 
°Ca II H and K chromospheric flux at 1 AU l|Browning et al.ll201fjh . 

Ca II H and K chromospheric flux at 1 AU l|Noves et al.lll984h . 
s Ca II H and K chromospheric flux at 1 AU jBaliunas et alj|l995ft . 

Ca II H and K chromospheric flux at 1 AU l lHenrv et al.lll99d) . 
'Ca II H and K chromospheric flux at 1 AU llHall et al.ll2007l , |2009 | ; lLockwood et alj|2007l) . 
J Ca II H and K chromospheric flux at 1 AU llCincunegui et al.ll2007T) . 
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Table 4. Iron Depletions and Differences between R(line) and Least-Squares Fits 



Star Spectral Type [Fe/H] Difference in dex of R(line) relative to Fits 



C IV C II 1 Mg II 



HR 4657 F5V -0.70 0.08 0.07 -0.12 -0.02 

X Her F8 V -0.50 0.28 0.28 -0.02 0.01 

r Cet G8V -0.43 0.15 0.18 -0.01 0.16 

36 Oph A Kl V -0.39 0.25 0.24 0.14 

Speedy Mic K2 V -1.49 0.44 

61 Cyg A K5 V -0.35 0.24 0.25 0.02 



Table 5. Least-Squares Fits to Lyman-a; Flux vs. Effective Temperature and Rotation 

Rate 



Distance 


-frot 


A 


B 


Mean 


RMS 


from Star 


(days) 






Dispersion(%) 


Dispersion(%) 


1 AU 


3-10 


0.37688 


0.0002061 


41.0 


73.6 


1 AU 


10-25 


0.48243 


0.0001632 


32.3 


42.1 


1 AU 


> 25 


-1.5963 


0.0004732 


85.0 


99.8 


HZ 


3-10 


3.9358 


-0.0004054 


34.8 


46.7 


HZ 


10-25 


4.5460 


-0.0005631 


37.4 


44.0 


HZ 


> 25 


3.5737 


-0.0004686 


43.1 


51.7 
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Fig. 1.— Plot of the ratio of the Lyman-a to C IV 1548+1550 A line flux vs. the C IV line 
flux at 1 AU. Included are stars between spectral types F5 V and M5 V, divided into three 
spectral type bins, the quiet Sun and the active Sun at two different times. The solar data 
are indicated by Sun symbols, and the L symbol refers to a star with low metal abundance 
[Fe/H] < —0.30. The solid lines are least-squares fits for each spectral type bin excluding 
the L stars and the Sun. T he ratio for a Cen A is closest to the solar ratios. The errors bars 
are 20% for stars using the IWood et al.l (120051 ) correction for missing Lyman-a flux or 30% 
for stars using the France et al. I J2012ah correction for missing Lyman-a flux. 
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Fig. 2. — Same as Figure 1 except for the Lyman-o; to C II 1335 A line flux ratio vs. the 
C II line flux at 1 AU. 
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Fig. 3. — Same as Figure 1 except for the Lyman-a; to I 1304.86 + 1336.03 A line flux 
ratio vs. the O 1 line flux at 1 AU. The small dispersion of the line flux ratios for the M 
stars may be unrepresentative given the small number of data points. 
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Fig. 4. — Same as Figure 1 except for the Lyman-a to Mg II 2796 + 2803 A line flux ratio 
vs. the Mg II line flux at 1 AU. 
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Fig. 5.— Plot of the Lyman-a to Mg II 2796 + 2803 A line flux ratio vs. the Mg II line flux 
at 1 AU for the MO to M5 stars. The solid line is th e least-squares fit to the data. The errors 



bars are 20% for stars usin g the 



Wood et al 



rt2005l ) correction for the missing Lyman-a flux 



or 30% for stars using the France et al.l ( I2012al ) correction for the missing Lyman-a flux. 
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Fig. 6. — Plot of the Lyman-a to Ca II 3933 + 3968 A line flux ratio vs. the Ca II line 
flux at 1 AU for the MO to M5 stars. The solar data are indicated by the Sun symbol, and 
the L symbol refers to a star with low metal abundance [Fe/H] < —0.30. The solid lines are 
least-squares fits for each spectral type bin excluding the L stars and the Sun. 
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Fig. 7. — Plot of the ratio of Lyman-a to X-ray flux vs. the Lyman-o; line flux at 1 AU for 
all stars in our sample with X-ray data. The solid lines are least-squares fits to the data for 
the F5 V-G9 V, KO V-K5 V, and MO V-M5 V stars. Data for the quiet Sun, moderately 
active Sun, and active Sun are indicated by the dotted circle symbols (from left to right 
respectively). M dwarfs stars are identified by name. 



-32 - 



3 



< 



X 

~1 



O 



o 



Circles: Exoplanet host stars 



2 







P(Rotation) 

— 3 days )K 
3-1 days O 
10 — 25 days 
P>25 days X 
P unknown + 




GO 



G5 



K0 K2 



K5 



10 



13 M5I 



7 000 



6000 5000 4000 

Stellar Effective Temperature (K) 



3000 



Fig. 8. — Lyman-a flux at 1 AU vs. stellar effective temperature. The stars are grouped 
according to stellar rotation period: ultrafast rotators (P rot < 3 days), fast rotators (3-10 
days), moderate rotators (10-25 days), and slow rotators (> 25 days). Rotation period is 
a rough measure of the magnetic heating rate in the star's chromosphere and corona. Host 
stars of exoplanets are circled and the quiet Sun is marked as a circled dot. Least-squares 
fit lines are shown for the fast, moderate, and slow rotators. 
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Fig. 9. — Same as Figure 8, except that the Lyman-a flux is evaluated in the habitable zone 
for an exoplanet. 



